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Abstract 

We study CP-violation in resonant three-body DPP decays, where 

PP stands for either Dtt, DKg, tttt or ttKs- Analogously to the — > Svr 
channel and the extraction of 2a, the first three channels are shown to mea- 
sure cos 2/3 in addition to sin 2/3, thus allowing to resolve the /3 — > ^ — /3 
ambiguity, while the DttKs final state leads to a measurement of 2/3 -|- 7. 
The B'^{t) D^D~TT^ channel via the interference between D** orbitally ex- 
cited resonances is taken as an example, although this Cabibbo-suppressed 
decay suffers from irreducible penguin uncertainties. Then two penguin-free 
and Cabibbo-dominant modes are proposed: B^^{t) D+D-Rs with the D** 
resonances, and B%t) D^tt+tt- with the D** plus the p. Finally, the 
D^TT^Ks channel with the Df and K* resonances provides a new 
clean method to measure the unusual angle 2(3 + ^. We present in all cases a 
crude estimate of the number of cos 2(3 (respectively 2(3 + 7) sensitive events. 
We show that this number is an increasing function of the resonance mass, a 
favorable situation compared to the more extensively studied three-pion Dalitz 
plot. However, the poor detection efficiency of the D mesons could pose a 
problem. As an annex and speculative application of these Dalitz plot based 
methods, the penguin-dominated B^{t) — > KsKsKi decay also measures 2(3. 
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Introduction 



There is little doubt that the CP-asymmetry in B ^ J/ipKs will be measured 
in the next years with an increasing accuracy and will provide sin 2/3 without 
sizeable hadronic uncertainties. However, measuring only sin 2/3 leaves ambiguities, 
in particular the ambiguity /? — *• | — /3. It might be argued that this is not a real 
problem since the overall constraint on the Unitarity Triangle (UT) already leave no 
room for this ambiguity, indeed we know that cos 2/3 > However, if [3 turns out 
to be not too far from 7r/4 an information on the sign of cos 2/3 will be important, 
and furthermore, the UT constraints only hold within the Standard Model, which 
obviously has to be tested with care in 5-factories. 

The analogous ambiguity a — > | — a was argued in to be solvable via the use of 
the time-dependent B^{t) —>■ tt+tt^tt'' Dalitz plot in the p resonance region. The in- 
terference between, say, the p+vr^ tt+tt^tt^ and the B^ p^7r+ ti^ti^ti^ 
decays and the use of the Breit-Wigner resonant phase allows to measure a term 
proportional to cos 2a and thus resolve the ambiguity. The partial understanding 
of the strong phase in the resonant region and the use of the isospin analysis fur- 
thermore provides some hope of separating the penguin contribution (which is far 
from negligible in this case and does introduce a significant shift on a) from the tree 
one 

One may wonder if the same technique may be applied to the measure of the 
angle (3. The straightforward generalization of the pvr — >■ Stt analysis would be a 
D*D — > D^D^TT^ one. However, the interest of the Dalitz plot analysis relies on the 
interference region, say, between the p"*" and p^ bands. This interference region is 
rather large thanks to the sizeable width of the p resonance. The D* is much too 
narrow (F < 0.1 MeV) to leave any hope. Therefore we considered an analysis of 
5° ^ D+D**- D+D-n^ interfering with 5° ^ D**+D- D+D-ir^ and with 
the CP-conjugate channels, where D** represents the orbitally excited positive parity 
charmed mesons (Z^g, Dl, Di, D2). In the final three-body state we can at will replace 
D hj D*. In this letter we present a preliminary study of this proposal. We find that 
in principle cos 2(3 can be measured this way, up to irreducible penguin uncertainties 
that cannot be removed via the Dalitz plot, due to the AJ = 1/2 isospin structure of 
the decay. 

We therefore present also two other proposals, free of penguin uncertainty and 
moreover Cabibbo dominant: B^ D^D**^ D^D^Ks interfering with B^ 
D+D**- D+D^Ks- This decay is the "class-I" analogous of 5 ^ J/ipKs. An- 
other interesting channel is B^ D**~it~^ D^pH^ir' interfering with B^ — >■ 
P'**+7r- ^ D^pTT+TT- and D%p^ D^pTr+TT", where D% means a neu- 

tral D decaying into a CP-eigenstate. For example, D^{D^) tttt, D^{D^) KK, 
D^{D^) KsiT^ etc. 

All above mentioned decay channels measure 2(3 since, in Wolfenstein parametri- 
zation [Q, the decay amplitudes are CP-invariant and the B^ — B^ mixing forwards 

^Throughout this paper, B means a Bd meson. 
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a exp[-2ip] phase. The Dahtz plots for B nD** + DK* D^n^Ks, on the 
contrary, contain the phase 7 of Vub and give a measure of 2(3 + 7. Such a measure 
would be very important, not only as an attempt at the measure of 7 but also to 
reduce further the angle ambiguities. Other methods to lift discrete ambiguities in 
CKM angles were discussed in [|, ^ . 

Unhappily it is not clear if the statistics will allow these analyses to be done in 
e^e" i?-factories (BaBar, Belle), mainly because of the small detection efficiency of 
the D mesons in the final states. However, except B —>■ D^D~7r^, all the presented 
B — s> DPP decays can be detected with charged particles only, which is a good 
starting point for hadronic machines (LHC-B, BTeV). 

This letter is organized as follows: in Section |I| we discuss the case of 5 — > 
D^D~TT^ explicitly and define an Effective Branching Ratio, which describes roughly 
the number of events generated by the Breit-Wigner interference effects and sensitive 
to the CKM angles. In Sections |, |, | we describe briefly the B D+D-Rs, 
B D^pTT^Ti^ and B D^n^Ks channels. In Section |^ a speculative idea on B ^ 
KsKgKi decay is presented. Finally in Section ^ we make some model calculation 
of the various effects described here. 



1 The B D+D-7T^ Dalitz plot 

1.1 The cos 2/5 dependence 

We will take the example oi B D^D^n^ just to be specific, and we neglect pen- 
guins. In Wolfenstein parametrization the weak decay goes through b — >■ ccd, and 
contains no CP-odd phase. The only phase in the problem will be the 2(3 phase of 
B^ — B^ mixing. Up to trivial angular rotations, two independent variables label the 
final state, which may usefully be taken as the Dalitz plot variables: 

where the relation 

s+ + s" + s° = m| + 2ml + (2) 

tells that there are only two independent variables, say and s~ . CP-eigenstates 
are on the line = s~ . Let us define 

A{s+,s-)=A{B''^D+D-7r^), s") = ^ Z^+Z^-yr"), (3) 

where B^ and B^ represent unmixed neutral B. The amplitudes in contain no 
weak phase, as already stated, they do contain unknown strong (CP-even) phases, 
which obviously depend on and s~ and are different in A and A except when 
s"*" = s~. 

The time- dependent amplitudes for an oscillating state B^(t) which has been 
tagged as a B^ meson at time t = is given by 

A{s+, s--t)= A{s+, S-) cos (^^) + te-^'^Ais^, s') sin 1^^) , (4) 
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and the time- dependent amplitude squared is: 



\A{s 




with 



Ge(s+,S-) 

Gs(s+,s-) 




(6) 
(7) 



(8) 



The transformation defining the CP-conjugate channel B^[t) —>■ D D^n^ is s+ s , 
A and 13^ -13. Then: 



Note that for simplicity the e and constant phase space factors have been omitted 
in Eqs. (|) and (|). 

Thus, if Gc(s"'',s~) is not antisymmetric in (s"^,s~) there is direct CP-violation 
(proportional to cos Am t); similarly, if Gs(s"'",s~) is not antisymmetric there is 
mixing- induced CP- violation (proportional to sin Am t). Neglecting penguins, no 
direct CP occurs, as A{s'^,s~) = A{s~,s^). Furthermore, if the final state was a 
CP-eigenstate, Gs would only contain the mixing-induced CP- violating sin 2/? term 
in (§). The CP-conserving cos 2/5 term is the new one. Measuring this term resolves 
the /3 — >■ I — /3 ambiguity But such a measurement necessitates some knowledge of 
the CP-even phase of A{s^ , s~)A*{s~^ , s~). In general the strong phases are utterly 
unknown a fortiori in a three-body decay. However the trick, first proposed in the 
three-pion case , is to assume resonance dominance, at least in chosen parts of the 
Dalitz plot, and to use the Breit-Wigner phase. 

We will consider the orbitally excited D** and D** resonances. The obvious 
problem here is that besides the Breit-Wigner phase, many other sources of strong 
phases can be thought of, for example due to the final state D** interaction with 
the D, etc. Nobody knows how to compute these phases. What we may however 
assume is that these background phases vary slowly under the resonances, where the 
Breit-Wigner phase varies quickly. Thus a fit with a constant background phase plus 
the Breit-Wigner phases should do the job. The second objection stems from the 
existence of many other Dtc {Dtt) resonances: higher excited charmed mesons. One 
cannot use too many parameters in the fit. It is assumed that one can minimize the 
effect of the unwanted resonances by selecting properly the domain in the Dalitz plot 
under the first orbitally excited D** and D** resonances. In our analysis we did not 
either consider, for simplicity, the charmonia DD resonances. A recent application of 

^It does not resolve the /3 — > tt + /3 ambiguity which requires some theoretical input 0. 
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these ideas in the non-CP studies of the D and Ds three-body decays may be found 
in §. 

In the following we will use for the masses (in GeV): 

ruB = 5.279, mo = 1.865, mo* = 2.360, mo* = 2.459. (10) 

Since the J = 1 orbitally excited mesons cannot decay into Dn (parity plus angular 
momentum conservation), we restrict ourselves to the J = 0, 2 ones. The J = 2, 
1)2(2459) has been observed and is narrow (23 ± 5 MeV) as expected since it decays 
through a D-wave. In our calculations we have taken r2 = 20 MeV. The Dq decays 
through an S-wave and is expected to be broad, which is enough to explain that it 
has not yet been observed. We will assume Fq = 150 MeV [||. 

For the sake of simplicity let us first consider only the contributions of the J = 
0, Dq and Dq resonances 



/(s+) = Breit-Wigner(D*+ ^ D+7r°) = — — — ^ (11 



s+ - m|,, + irriD'To 

A{B^ D+Dl- D+D-7i°) = Ti /(s") (12) 

A{B^ D*+D- D+D-7i°) = T2 f{s+) (13) 

A{B^ D-DI+ D+D-7i°) = Ti f{s+) (14) 

A{B° D*-D+ D+D-7i°) = T2 f{s-) (15) 

where BRq is the branching ratio Dq'^ — >■ D~^7t^, is the momentum of the D in 
the Dq rest frame, and where we have used the CP-identities in Wolfenstein phase 
convention 0: 

M{B^ D^D*-) = M{B^ D-D*+) = Ti (16) 
M(5° ^ DI+D-) = M{B^ D*-D+) = T2. (17) 

Our convention is to write first, in the final state, the meson which is directly coupled 
to the initial meson via a current form factor, and second, the one which is emitted 
by the W meson. The numerators in ([Tl|) have been fixed by imposing that in the 



narrow width limit the integrated three-body decay width for B^ D~^Dq~ — >■ 
D^D^Ti^ coincides with the two-body one B^ D^Dq^ times BRq. For resonances 
of spin J there is an additional multiplicative factor of Yj {9* , 0) / \^ , where 9*,(f* 
are the spherical angles of the decay products in the resonance rest frame. Angular 
momentum conservation selects m = and hence there is no dependence on ip*. 
Defining 

Gtot = IT1P + IT2P (18) 



'''It should be stressed that in a real analysis, one would take into account as many resonances as 
possible. See, for example, ref. |7|. 

^Any constant phase in the normalization of the Breit-Wigner functions can be incorporated in 
the definition of Ti and T2. 
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R = (\T,\'-\T2n (19) 

^tot 



D = = -^\Ti\\T2\ (20) 

^tot 

S = Arg(riT2*) (21) 



we obtain: 



Go(s+,s-) = G,ot{\f{s^)\' + \f{s')\' + 2D cos5Re(f{s+)f{s'r)} (22) 

G,{s+,s-) = Gtot{i?[|/(5+)p-|/(s-)P] -2Dsin5Im(/(.+)/(s-)*)} (23) 

G,{s+,s-) = Gtot{-/^[sin(2/5 + 5)|/(s+)|2 + sm(2/3-5)|/(s-)|2^ 

- 2sm2/3Re(/(s+)/(s-)*)-2i?cos2/5Im(/(s+)/(s-)*)} (24) 

The important point is that by measuring the kinematical and time-dependence 
of the three-body decay, one extracts from (p^-p^) all the two-body complex ampli- 
tudes (^^-|T^), as well as 2/3, without discrete ambiguities in the general case. This 
is true whatever the number of contributing resonances is, as long as we are able to 
deal with the multi-parameters fit 0. 

The amplitudes Ti and T2 are not in general related. From (0), R is thus typically 
of order 1. Then one is lead to use the last term of the rhs of (|2^) to measure cos 2f3. If 
by some symmetry or by chance R is small, D is of order 1 (in fact both R and D are 
often of order 1). The first two terms in the rhs of (|2^) may then be used. There the 
cos 2(3 term comes also from an interference between strong phases, but instead of the 
Breit-Wigner phases, it uses the strong phase S of the decay B D**D. This phase, 
which vanishes in the factorization approximation, might be small ~ l/N^. To get it 
one has to use eqs. (^2]) and ( p3D where it is measured thanks to the Breit-Wigner 
interference. This shows that the determination of cos 2/3 is overconstrained. 

To summarize, the extraction of both the sine and cosine of the CKM phase rely 
on the important hypothesis (hopefully very reasonable) stated as follows: 

Over the Dalitz plot, no other CP -even phases than a la Breit-Wigner, 

or constants, are allowed. 

In the B D^D^tt^ channel, R is indeed presumably close to 1, as will be 
discussed in a model in Section ||, because the B ^ D transition is strongly favored 
over B — > , while the current matrix elements to the vacuum are presumably of the 
same order of magnitude. On the contrary, for the D2 contribution, one has R = —1 
exactly in the factorization approximation as the D2 does not couple to the current. 

From (|2^), one sees also that whatever the values of R and D are, a sin 2/? term 
is present. This term will not compete with the measurement of the B J/ipKs 



^This statement is valid assuming the resonances are well known (mass, width, spin...). This 
could be a difBcult experimental task, particularly for the Dq. 
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asymmetry, but it can be used either as a cross-check or as an additional input which 
would help a multi-parameters likelihood fit. 

Now, the factor f{s~^)f*{s~) will be studied in some detail in the coming sub- 
section. It will be seen that the wanted effect grows like the resonance width and 
like its mass squared. This is why rather larger and heavier resonances are better. 
They should however not be too large in order to emerge from the background noise. 
A priori, the non-yet observed Z)q resonance is large enough ~ 150 MeV, while the 
known D2 is rather narrow ~ 20 MeV. 

In this section we have neglected the penguin diagrams, although this is not really 
legitimate. Unhappily, the Dalitz plot analysis does not allow to separate penguin 
contributions from the tree ones, as it can be done in the Sir case because for 
the b — > ccd transition , trees and and penguins have the same A/ = 1/2 isospin 



structure |10| . Of course it might be argued that since the aim is to resolve a discrete 
ambiguity on (3, a small penguin induced error is not a problem. But a penguin free 
channel would be safer, and this is the theme of the next sections. Before that let us 
in the next subsection work out a crude estimator of the number of cos 2(3 sensitive 
events in the Dalitz plot. 

1.2 Effective Branching Ratio 

Throughout this letter we are interested in Breit-Wigner interferences, either to gen- 
erate a cos 2/? term, as in DDn, DDK and Dun final states (see (p^), or to produce 
2(3 + ■J dependence, as in DnK final state where all the CKM angle dependence is 
proportional to the Breit-Wigner interferences (see (^41) ). 

Let us call fi the Breit-Wigner of the resonance Ri the mass, width and spin 
of which is respectively m^, Fj, Jj. Recall eq. (|2^): the interference term fiYif2Y2 
describes the fluctuations of the number of events which are due to cos 2(3. What is 
the size of this term 7 Yi = \/AtiYj^{0*, 0) with 9* in the rest frame of resonance Ri 
corresponds to the angle between the resonance decay product momentum and the 
momentum of the B. When the resonance is in the channel, cos^* depends linearly 
on s~, and vice- versa. In practice the resonance crossing region corresponds to 9* ^ 
which tends to enhance the interference effect for large spin of the resonances. To 
get a flavor of the result, let us first consider the simpler case of Ji = J2 = 0: 



iVi2 = / ds+ds-\\m{h{s+)f;{s-))\ = STrmimsA/BRiBRs^^ x (25) 
J V P1P2 

where 



BRj is the branching ratio of Ri in the considered decay products, and where p* is 
the final momentum of these decay products in the resonance rest frame {p* ranges 
from 300 to 500 MeV in the considered cases). We end up with a term proportional 



to S'Kmim2\ TiV2/{piP*2) times a dimensionless integral. The latter depends on the 
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masses and widths through the integration bounds which are typically of the form 
{smax — ''TT''i)/{^ifni)- A rapid inspection shows that the dimensionless integrals cannot 
depend on the phase space parameters and widths more than logarithmically. So the 
main dependence is included in the dimensionfull factor: 



Fi2 = 87rmim2 JbRiBR2^4^. (27) 
V P1P2 

We immediately see that heavy resonances will be favored. This can be partially 
understood as the width of the resonance band in the Dalitz plot is mT and not 
r simply. The results of exact numerical integration of (pSf ) will be given for each 
channel in Tables |1], 0, ^, ^. It happens that the dimensionless integral is larger when 
the resonance crossing is inside the physical domain than when it is not. In Figure |I|, 
we have plotted the shape of two representative phase spaces, namely B — > DDn and 
B —>■ Dun. In the former case the j cross well inside the physical domain, 
and I12 is typically ~ C*^* x log [{smax — ''^1) / ij^i^ i)] ~ 5 — 20, while in the latter the 
£)** I jj** cross well outside the physical domain, and /12 is of order 1. This is expected 
since we are looking for an interference effect between these two resonances, but it 
is noticeable that even in the latter case, there is a sizeable dependence on cos 2/? 
located far from the resonance interference region. This feature, already seen in the 
three-pion case 0, |^, is due to the fact that the Breit-Wigner formulae, when they 
are not squared, as in (^5]), decrease slowly away from the resonance center. However, 



as already discussed, the domains too far away from the resonance regions need much 
care since other resonances, not considered in the fit, may be dominant in the latter 
regions. When the resonances do cross in the physical domain, a very conservative 
estimate of the cos 2/5 sensitive events is equivalent to take J12 = 1 in (pSf), or 



/i2,cons ^ 47ry'F0 (0, 0) (0, 0) = ^J{2.h + 1)(2J2 + 1) (28) 

in the case of resonances of spins Ji, J2. When they do not cross, the conservative 
estimate is simply 0. 

Now we are able to define an Effective Branching Ratio (BRegf) for the contribution 
of Dl to the cos 2/5 term (see eqs. (^) and (^Sf )): 



BR, 



cff 



BR(5° ^ DI+D-) + BR(fi° ^ D+DI-] x Rx , (29) 

J 32n^mB < p > 



with the two-body BR calculated in the zero width limit (in the following, we use the 
factorization assumption to compute these BR) and 



Pi 




'+P2IT2 


2 




Ti 


2 + 


T2 


2 



<p>= ,^ ,„ ,^ . (30) 



Pi (resp. P2) is the momentum in the rest frame of the B of the resonance 1 (resp. 
resonance 2). 
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Figure 1: To the left, the shape of the B DDn phase space (the straight lines 
represent the position of the D** resonances). To the right, the same for B Dtttt 
(the third straight line represents the p). It is seen that the jjj** resonances cross 
inside the DDtt physical domain, and outside the Dtttt one. 





B\t) D+D 




= 2/3 


Crossing 


BReff for cos$ 


Total BR 


Detection Efficiency 




(0.2-3) X 10-5 


21x10-5 


0.005-0.01 


D*+/Dl- 


(0.01-0.04) X 10-5 


1.9x10-5 




Total 


(0.1-3) X 10-^ 







Table 1: Effective Branching Ratio for cos 2/3 in the B — > D^D 7r° channel. 



The fact that oscillations generate the cos 2/3 factor is obvious from (|^) since 
Gs multiplies sin Ami. We must expect an additional suppression due to these 
oscillations. However, to get an idea of this suppression one integrates over time 
exp(— Ft) sin Ami, which gives x/{l + x^) ~ 0.48 where x = Am/F ~ 0.73. We will 
neglect this factor 1/2 in view of the crudeness of our estimate. 

For the B D^D^n^ Dalitz plot, the Effective Branching Ratios concerning the 
contributions of Dq and D2 are shown in Table |l|. The lower bound corresponds to 
the conservative estimate for J12 (see (|5)) while the upper bound corresponds to 
the full numerical integration of (^). The two-body non-leptonic branching ratios 
have been calculated in the factorization assumption, as explained in section ^ The 
last line of the table sums all the contributions taking into account the detection 
efficiency of the particles in the final state. 

For comparison, the Effective Branching Ratio for the measurement of cos2q; in 
the 371 case is about 10-^. In this case the / p~ crossing lies just at the border of 
the physical domain. From Table we see that the possibility of measuring cos 2(3 
m B ^ D^D^Ti^ channel relies strongly on the contribution of the Dq (because the 
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B%t) D+D-Ks ^ = 2(3 


Crossing 


BRoff for cos $ 
(0.6-6) X 10"^ 


Total BR 
35x10"^ 


Detection Efficiency 
0.005-0.01 


Total 


(0.3-6) X 10-6 





Table 2: Effective Branching Ratio for cos 2/? in the B D^D Ks channel. 



transition B^ D^Dq^ dominates over B^ D^D2 and B^ D2^D^) and thus 
on our ability to describe accurately this resonance. In conclusion, this channel seems 
to be not accessible to an e^e' i?-factory and the presence of a 7r° is problematic 
for an hadronic machine. However, it is straightforward to generalize this idea to 
B — > D'^*'>^ D^*^~ Ti^ channels, for which the detection efficiency will be much more 
favorable. 

In the following, we will consider cases where the two crossing resonances are not 
necessarily CP-conjugate. The calculation of BRefj follows the same line and only the 
numerical results will be presented in the next sections. 



2 B ^ D+D-Ks 

At the quark level, the transition generating this channel is 6 — ccs, that is, the 
same as in B ^ J/'^Ks- This has two important consequences: First, it is a class-I, 
Cabibbo-favored channel (proportional to |aiVcfeVcs|). Second, the penguin pollution 
is suppressed [] by |\46V^s|/|KfeKs| ~ 2 x 10~^. The kinematics of this decay is 
very similar to i? D^D~7r^, replacing the D** resonances by the D**. Due to 
its flavour structure, the D** contribute only when connected to the W propagator. 
Thus the contribution of the D*2 is expected to be very small (it is zero in the 
factorization approximation), and the non-yet observed D*q should dominate the 
decay. We assume mD*^ = 2.5 GeV; for the width, we take Tgo = 150 MeV and assume 
that DK saturates the decay p. The equations (|2^^^ apply, with R = 1, D = 0. 



The calculation of the Effective Branching Ratio, summarized in Table 0, is very 
encouraging (~ 10"*^), although it has once more to be corrected by the D meson 
detection efficiency. In any case, this channel should be about 20 times more favorable 
than the B D+D-JC^ one. 

Finally, this analysis is not expected to be generalizable to the analogous decay 
with one (or two) D* in the final state, as the D*q does not decay in D*K, the D*^ 
is expected to be below the D*K threshold, the D^i is too narrow ( F < 2.3 MeV), 
and the branching ratio for the D*2 is expected to be suppressed. 



^We have used CKM unitarity to incorporate the |Vt6l4s| contribution in the "tree" part propor- 
tional to |Vc&Vcs| and in the "penguin" part proportional to |KibKis|- 
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3 B D^pTT+TT- 

By D(^p we mean a meson which eventually decays into a CP-eigenstate, for 
example Kst^^-, vr+vr^, vr^vr^, K^K~ K^K^, i.e. a few % of the are used. How- 
ever, compared to the situation in the preceding sections, this loss in statistics is 
similar to the two D mesons detection efficiency. Furthermore, as -B — > D~^D~Ks, 
B — *• D^pTT+TT" is Cabibbo-favored (proportional to The very noticeable 

advantage of B D^piT^Tf' is that there is no penguin contribution at all. Indeed 
the weak decay at the quark level is of the type b cud and obviously no penguin 
operator can contribute to such a transition with four changes in flavor. 

The analysis proceeds then in a way parallel to the one in the preceding sections. 
The two-body channels will be: 

B°^D**+7i' and B^ ^ p^D^, (31) 

and their CP-conjugates. The channels P° tt^D**^ and P° p^D^ are doubly- 
Cabibbo-suppressed and we will forget them. Hence the ratio R (|19|) relevant for 
the crossing between B^ — > D**^7i^ and P° D**^7i^ is equal to 1. The branching 
ratios recently measured by CLEO for B Din and B D^^t^ are of order 
10-3. 

It should be noted here that the crossings are outside the physical do- 

main. We then better consider the D** / p crossings which are inside the physical 
domain. The problem now is that it is difficult to estimate accurately the ratio: 

\A{B'^p'D')\ 
\A{B^ ^ D**+TT-)\' ^ ' 

A{B^ D**~^7i~) is suppressed by the ri/2 or T3/2 form factors, while ^4(5° p^D^) 
is color-suppressed. The ratio (|32D is expected to be of order 1, as happens in the 
factorization assumption, with the form factors computed in our model described in 
section ||. 

The Table |^ shows the Effective Branching Ratio for this channel, which is compa- 
rable to the one corresponding to the measurement of cos 2a through the Stt channel. 
In spite of the relative narrowness of the D^, the -D2/P crossing is predicted to be one 
dominant contribution to the measurement of cos 2(3, thanks to the spin effects, the 
width of the p and the large branching ratio for B ^ D'^t^ compared to P — > DqT[. 
This is nice, because the D2, which is already seen in semi-leptonic and non-leptonic 
B decays, will be well known in a close future. 



^Proportional to iKibVcdl, they would introduce a very small dependence on the angle 2/3 + 7, 
and could in principle be incorporated in the analysis. However, the B D^n^ Ks channel is much 
more favorable in this respect. 
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Crossing 


BRcff for cos $ 


Total BR 


Detect fon Efficiency 




(0.1-0.4) xlO^^ 


16x10^^ 


0.005-0.01 


d;+/d;- 


~ 0.4x10^^ 


120x10^^ 






(0.1-2) X 10^5 


32x10-5 






(0.2-4) X 10^5 


130x10^5 




Total 


(0.05-1) xlO-t* 







Table 3: Effective Branching Ratio for cos 2/5 in the B D^p-n^n channel. 



4 B ^ D^TT^Ks 

In the case of the D^n^Ks final state, the intermediate resonant channels are 

B^ D+K*' ^ B\ B^ ^ Ti-D*;+ ^B\ B^ ^ KsD**^ ^ 5°. (33) 

The K* contributes only to diagrams proportional to | VusKfe] (= 0{X^) in Wolfenstein 
parametrization) while the D** contribute only to diagrams proportional to |KsV^b| 
(= 0{X^) too). Since, at the quark level, the weak decays are b — * csu and b —>■ use, 
there is no penguin contribution in this channel with four fiavor changes. This very 
interesting channel measures cos(2/5 + 7) and sin(2/5 + 7). Indeed, in Wolfenstein 
parametrization, the CKM matrix element Vub contains the phase 7 while the B^ — B^ 
mixing is proportional to exp(— 2i/3). To our knowledge, the only other ways to 
measure sin(2/9 + 7) is to look at the very small CP-asymmetry in the dominant 
channel B — > D{D)7t, the larger CP- asymmetry in the color-suppressed channel 
B KsD^{D^), or in some semi- inclusive b — > ucs + cus decays |Tl|; no other 



method to measure directly cos(2/5 -|- 7) has been proposed in the literature. 

As this channel is a little different from the others presented in the previous 
sections, we rewrite equations (|22|-p4D for the D~^n~Ks final state. Neglecting for 



simplicity the neutral mode KsD**^, the following channels contribute: 

Anr = A{B' ^ n- DU) = Tn^^/^ (34) 

Ak* = A{B^ ^ D+K*') = Tk* (35) 

where we made use of the fact that the tensor (J=2) D*2 particle cannot couple to 
currents from the vacuum, and thus A{B^ 71^0*2) is suppressed. Defining: 

Gtot = \TD:f + \TKA', (36) 

\Tk'\'), (37) 



tot = 




|2 1 \rp |2 

1 + l-LR*] , 


R = 


1 

Gtot 


{\To..f - 


D = 


Vi- 


Gtot 


5 = 


Arg( 





\Td:,\\TkA. (38) 

(39) 
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B%t) - 


D^n^Ks $ = 2/? + 7 




Crossing 


BReff for cos $ 


BReff for sin $ 


Total BR 


Detection Efficiency 


d:^/k*- 


(0.04-0.2) X 10"^ 


(0.04-0.4) X 10"^ 


27 10"^ 


0.05-0.1 


dU/k*o- 


(0.1-2) xlO^^ 


(0.1-2) X 10"^ 


27 10"^ 




Total 


0.07-2 10"" 


0.07-3 10-" 







Table 4: Effective Branching Ratio for cos(2/3 + 7) and sin(2/3 
D^-k^Kq channel. 



7) in the B 



we get 
Go 

Ge 

G« 



G 



tot 



(7+ = Breit-Wigner(D*(]" 
/i_ = Breit-Wigner(i^* 



l + R, ,2 1-i?,, ,2' 



D+Ks) 
KsTi-). 



\9+r + 



Gt, 



l + R. ,2 |2 



(40) 
(41) 



(42) 
(43) 



Gtnt-D 



cos((5 + 2f3 + 7)Im - sin((5 + 2/3 + 7)Re • (44) 



The B^{t) D^Ti^Ks time-dependent decay allows a measurement of the angle 
5 + 2/9 + 7, while B^(t) D^ti^Ks leads to 5 — 2/3 — 7. In the end, one obtains 
separately the CP-violating phase 2/3 + 7 = 7r + /3 — a and the CP-conserving phase 
5. Once 2(3 is known (cf. above), this could be a very interesting measurement of 7. 

The contribution of the i^'*(892) resonance to the 2/3+7 terms suffers from its rela- 
tively small mass and width and more importantly from the fact that the K* (892) / D*q 
crossing is outside the physical domain. Thus it is important to take into account 
higher K* resonances such as the broad (287 MeV) scalar irQ(1430) (the contribu- 
tion of the although interesting, vanishes in the factorization assumption): these 
crossings lie just at the border of the physical domain. From Table ^ the mea- 
surement of both sin(2/3 + 7) and cos(2/3 + 7) from DtiK analysis should be a task 
comparable to the measurement of cos2q; from 37r channel, although some care has 
to be taken deciding which resonances should be cut or retained. 



5 Another cos 2/3 channel: B KsKsKi 

Finally we present as a speculative idea the B — ^ KsKsK^ Dalitz plot where the 
resonant channels to be considered could be: B — (pKs and B — > /o(980)-ft'L. From 
Bose statistics decays into K^Ks in P-wave, while /o(980) decays into KsKs or as 
well KlKi in S-wave. Only penguin diagrams contribute to these decays, which are 
real in Wolfenstein parametrization, except for suppressed long-distance penguins, 
and thus the mixing angle 2j3 is measured by the Dalitz plot. The narrowness of the 
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(j) resonance may look worrying for this to work, the nature of the /o(980) is unclear 
{qq or "molecule"), and anyhow a closer scrutiny of this idea is needed. The crude 
estimator described in section L2 gives in this case an Effective Branching Ratio of 
order 10~^ or less (assuming B cpKs ~ -B ^ IqKl ~ 10~^). 



6 A simple model for B D^D tt^ 

First we must repeat that our aim in this letter is mainly illustrative. We will use the 
factorization hypothesis for the non-leptonic decays. We will assume Heavy Quark 
Symmetry for the wanted heavy-to-heavy form factors and leptonic decay amplitudes. 
Finally we take the form factors as computed recently in a class of relativistic quark 



models ||T3|. Let us first give a brief account of the form factors and parameters we 



have used 



9. 



^H=0.9(-^)', ri/2H = 0.23 (-^)''\ r3/2H = 0.54(-^)', (45) 
\l + wj \l + wj \l + wj 



were we use the form factors in the notation of |jT^ and the indices 1/2 (3/2) refer to 
the angular momentum j of the light quanta in the D** mesons, which is conserved 
in the heavy mass limit. 

For the leptonic decay constants we take 

Id = fn,,, = 200 MeV, fo,,, = 0. (46) 

The vanishing of /D3/2 is a general consequence of the heavy mass limit |T5[ . 

The non-leptonic decays D**~^D~ and 5° —>■ D^D**~ are computed via the 

factorization hypothesis: 

M(E0^Z}+D*~) = Kd,a^f^,^'^\mB-mD)^^^^E^{w + l)i{w), (47) 
M(B0^D*+D-) = Kd,a^2fD{mB + mD*).JmB^{w-l)ry2{w), (48) 

M{B''^D*+D~) = Kai>aiV2fn{mB + mD*),/mB^ rriB {w^ - l)r3/2(w), (49) 
M(5° ^ D+D;-) = 0, (50) 

where K^;, = ^V*^Vcb and the last equality follows from the fact that the D2 has no 
matrix element to the vacuum due to /D3/2 = IH • 



The three-body amplitudes are obtained by multiplying the above amplitudes 
with the Breit-Wigner formula: 



mD*,/87rBRo^ 

A{B° ^ D+D*^- D+D-7r°) = M(5° ^ D+D*-^ " ^ 



s - mi* + immFo' 



^In (p5|), and contrary to ref. ||T^, we include a factor 0.9 in £,{w) to take into account perturbative 
and non-perturbative corrections to ^(1), thus avoiding overestimates of the rates. 
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with w = -'B + mn-s- 
2mBm£) 

mz)., /SttBRo^ 

^ ^ D+D-7r°) = M{B'' ^ D*+D-)— ^ 

s+ — mf)* + im D* To 

with w = ■■^Bj'^ - "^D ^^2) 



mg + s"^ — mj 
2mBmD* 



^0 



mD*./87rBR2^h(s 
yl(B° D;+D- D+D-n"^) = M(5° ^ D*+D-^ ' ^ 



with ^ = "^^^^ (53) 



m|j + s"*" — m 
2mBmD. 



2 



where (P2) is the momentum of the final mesons in the rest frame of the decaying 
Dq (-D2). BRo is the branching ratio Dq'^ D~^7t^, i.e. just a Clebsch-Gordan 
coefficient if we assume that Dir saturates the Dq decay: BRq = 1/3. BR2 takes into 
account an additional ratio of 70% of {D* DTr)/{D; Dn + D*Tr): BR2 = 0.23 
and finally the function 

= — (3cos2r-l) (54) 



2 V" 

is proportional to ¥2(6*, 0* = 0) function describing the D-wave decay of the D2, 9* 
being the angle of the D meson with respect to the B momentum in the D2 rest frame. 
The Breit-Wigner functions are normalized, as indicated in section |T|, to recover the 
two-body decay widths in the narrow width limit. Notice that in (pl|-^), the two- 
body amplitudes do depend on (s"*", s~) via the dependence of the form factors on w. 
It is often assumed in the fits ||2|, 0, ^ that the two-body amplitudes are constant, 
and admittedly it is difficult to avoid such an assumption in the present status of 
our knowledge, but it must be asserted that in reality they are only approximately 
constant, as we find in our model, as would be imposed by analyticity, etc. This is 
a source of theoretical error in the extraction of the CKM angles: its estimation is 
beyond the scope of this letter. 
Then 

A{s-,s+) = A{s^,s-) = A{B'> ^ D+D*^- D+D-n'') (55) 
+ A{B° -> D*+D- D+D-Ti^) + A{B^ DI+ D+D-ji^) 

and the differential branching ratio is given by 



dBR{B\t) ^ D+Z}-7r°) = l-^-}L_rB\A{s^,s--t)\^ds^ds- (56) 



where tb = 1.56 10^^^ s is the neutral B meson life time and \A{s~^ , s']t)\'^ is given 
by (P). The differential branching ratio dBK/ds^ds" as a function of the Dalitz 
point (s+, s~) is plotted on Figure 0. It is seen clearly that the cos 2/5 term produces 
a significant effect. 
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Figure 2: B^{t) D^D^n^ from our simple model. To the left, the differential 
branching ratio ^dj ) for Amt = 7r/2,/3 = 7r/6. To the right, also for Amt = tt/2 
the differential branching ratio j ^dj ) for j3 = tt/Q minus the one for (3 = 27r/6. This 
difference isolates the cos 2/5 term. 

Conclusion 

We have proposed several three-body Bd decay channels from which cos 2/5 can be 
extracted. When mediated by heavy resonances, the interference effects should be 
more important than in the light systems. 

The B D^D^Ks mode is the most promising as it is Cabibbo-dominant and 
almost penguin-free. For the same reasons B D^pn^Ti' is very interesting. Al- 
though not competitive in its naive version, the B D^D^tt^ channel should be 
directly generalizable to B ^ D*^D*~ti^ and thus would benefit from a bigger de- 
tection efficiency, though it suffers from irreducible penguin uncertainties. 

A new clean method to measure 2/? + 7 is described, which requires the analysis 
of the penguin-free B D^tt^Ks Dalitz plots. 

All the above charmed three-body decays might be detected at an e~^e~ 5-factory. 
The hadronic machines will be also a good place to look at the three most interesting 
channels, namely B D+D-Rs, B D^-k+ti- and B D^n^Ks- 

Finally, we notice that the penguin-induced B KsKsKl mode could also 
measure cos 2(3 although the branching ratio might be too small. 

As shown by the huge uncertainty of our crude calculations, further work is 
needed, along the line of Sir studies, to demonstrate the feasibility of these anal- 
yses. 
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